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Abstract 

The spherically symmetric dust model of LemaTtre-Tolman can describe wormholes, 
but the causal communication between the two asymptotic regions through the neck is 
even less than in the vacuum (Schwarzschild-Kruskal-Szekeres) case. We investigate the 
anisotropic generalisation of the wormhole topology in the Szekeres model. The function 
E{r,p,q) describes the deviation from spherical symmetry if dr-E ^ 0, but this requires 
the mass to be increasing with radius, drM > 0, i.e. non-zero density. We investigate 
the geometrical relations between the mass dipole and the locii of apparent horizon and 
of shell-crossings. We present the various conditions that ensure physically reasonable 
quasi-spherical models, including a regular origin, regular maxima and minima in the 
spatial sections, and the absence of shell-crossings. We show that physically reasonable 
values of dfE ^ cannot compensate for the effects of dr-M > in any direction, so 
that communication through the neck is still worse than in the vacuum. 

We also show that a handle topology cannot be created by identifying hypersufaces in 
the two asymptotic regions on either side of a wormhole, unless a surface layer is allowed 
at the junction. This impossibility includes the Schwarzschild-Kruskal-Szekeres case. 
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04.20. Gz, spacetime topology &i causal structure 
04.40.Nr, spacetimes with fluids or fields 
04.70. Bw, classical black holes 
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1 Introduction 



The Szekeres metric is a dust model, which has no Killing vectors but contains the 
LemaTtre-Tolman (LT) model as the spherically symmetric special case, which itself contains 
the Schwarzschild-Kruskal-Szekeres []TB|, |23|] manifold as the vacuum case. As with the LT 
model, it is written in synchronous coordinates, and the particles of dust are comoving. The 
constant time slices are foliated by 2-surfaces of constant coordinate r, which have 2-metrics 



of spheres, planes or pseudo-spheres, depending on the value of parameter e. See |]T7| for a 
review of its known properties. 

Despite the inhomogeneity of the model, and the lack of Killing vectors, any surface of 
constant coordinate 'radius' r in the e = +1 case can be matched onto a Schwarzschild vacuum 
metric [^, and any surface of constant time t is conformally flat [Q]. 

We here investigate the topological and causal properties of the quasi-spherical case, 
e = subject to the requirements for a physically reasonable model. Reasonability require- 
ments include, well behaved metric components, non-divergent density and curvature, regular 
spherical origins, regular maxima and minima in the spatial sections, and prohibition of shell 
crossings. Choosing well behaved coordinates also assists in avoiding the confusion of coordinate 
singularities. 

Studying such models of low symmetry is important, so that one can check which properties 
of spherically symmetric investigations of cosmology and gravitational collapse are preserved, 
and which are not. 

The subjects studied in this paper in some detail are the following: 

1. The dipole-like variation of mass-density; the locus of its poles and of the equator, and 
the images of the equator under the Riemann projection. 

2. Conditions for regularity of the geometry at the origin R = 0. 

3. Intersections of the shell crossings with the surfaces of constant (t, r), and conditions 
for avoidance of shell crossings. 

4. Conditions for regular maxima and minima at necks and bellies. 

5. Conditions for a handle topology of a t = const space, and the impossibility of preserving 
this topology during evolution of the model. 

6. Apparent horizons - their shape, intersections with the surfaces of constant (t, r), 
relations between these intersections and those of shell crossings, and with the dipole equator, 
location of an AH with respect to the R = 2M hypersurface, the intersection of an AH with a 
neck. 

7. The impossibility of sending a light ray through the neck so that it would emerge from 
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under the AH on the other side. 

8. Numerical examples of light paths traversing the neck and of those going in its vicinity. 



2 The Szekeres Metric 



The LT-type Szekeres metric isQ: 



where ' = d/dr, e = ±1, and / = /(r) > — e is an arbitrary function of r. 
The function E is given by 

E(r, p, q) = A{p^ + g2) + 25ip + 2^2? + C, (2) 

where functions A = A{r), Bi = Bi{r), B2 = B2{r), and C = C(r) satisfy the relation 

4{AC - Bl -Bl)=e , e = 0, ±1, (3) 

but are otherwise arbitrary. 

The function R = R{t,r) satisfies the Friedmann equation for dust 

R^-'-f^f, (4) 

where " = d/dt and M = M{r) is another arbitrary function of coordinate "radius", r. It 
follows that the acceleration of R is always negative 

■■ -M 

fl = ^. (5) 

Here M(r) plays the role of an effective gravitational mass for particles at comoving "radius" r. 
For e = +1, it is simply the total gravitational mass within the sphere of radius r. We assume 
M > and i? > 0. In /(r) represents twice the energy per unit mass of the particles in 
the shells of matter at constant r, but in the metric (|l]) it also determines the geometry of the 
spatial sections t =constant (c.f. [|10]). The evolution of R depends on the value of /; it can 
be: 

hyperbolic, / > 0: 

M 

R = y(coshr/-l), (6) 
f/'a{t - a) 

{smhr]-r]) = — , (7) 



■"■The results presented in Ref. contain a few misleading typos that were corrected in Ref. The 
notation used here does not follow the traditional one 
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parabolic, f — 0: 



2 

R = m|-, (8) 
rf ^ a{t - a) 

6 M ' ^ ^ 

I.e. R = V^j ' 



elliptic, / < 0: 



M 

R = ^(l-cos77), (11) 

(r^-smr^) = ^ ^, (12) 

where a = a(r) is the last arbitrary function, giving the local time of the big bang or crunch 
R = and a = ±1 permits time reversal. More correctly, the three types of evolution hold 
for f/M'^/^ >,=,< 0, since / = at a spherical type origin for all 3 evolution types. The 
behaviour of R{t, r) is identical to that in the LT model, and is unaffected by (p, q) variations. 

A more meaningful way to write E is 

.(...)^f{(^)%(^)%.), (13) 
where 5" = S{r), P — P{r), and Q — Q{r) are arbitrary functions, and 

^-2S' Z5 ' ^'-^ ' Ys ■ ^^^^ 



The metric component 
is actually the unit sphere, plane, pseudo-sphere in Riemann projection: 



(15) 



e = = cot f ^1 cos(0) , ii^ = cotf^'lsin(0), (16) 



(p-P) /2\ {q-Q) (2 



s \e 



(9 


-Q) 




s 




-Q) 




s 




-Q) 



2 



cos(0) , (17) 



6 = -1 ^^-^ = coth (^^ cos(0) , i^-^ = coth (^^) sin(0). (18) 

It seems reasonable to expect S > but it is not obviously impossible for 5" to reach or pass 
through zero. 

The factor e determines whether the p-q 2-surfaces are spherical (e = +1), pseudo-spherical 
(e = —1), or planar (e = 0). In other words, it determines how the constant r 2-surfaces foliate 
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Figure 1: The Riemann projection from {0,(p) to {p,q) coordinates for spheres & two- 
sheeted hyperboloids. The diagrams show only the = 0,7r section, i.e. the q — Q section. 



the 3-d spatial sections of constant t. The function E determines how the coordinates (p, g) 
map onto the unit 2-sphere (plane, pseudo-sphere) at each value of r. At each r these 2- 
surfaces are multiplied by the areal "radius" R = R{t,r) that evolves with time. Thus the 
r-p-q 3-surfaces are constructed out of a sequence of 2-dimensional spheres (pseudo-spheres, 
planes) that are not concentric, since the metric component depends on p and q as well as 
r and t. 

The (p, g)-coordinates in the cases e = +1 and e = have the range (— oo, +oo). In 
the case e = — 1, the parametrization ( |13D does not cover the subcases A = and C = 
(these subcases cannot occur with e > because of (|^)). Coming back to (^, we see that, 
for e = — 1 and A ^ 0, E \s zero when 

(p + B^/Af + iq + B,/Af = l/iAAY, 

E is positive for p and q outside this circle, and is negative for p and q inside it. Fig. |1| suggests 
that with e = — 1, we should rather take {—E) as the metric function so that p and q have finite 
rather than semi-infinite ranges. However, both the E > and E < regions are Szekeres 
spacetimes because they are mapped one onto another by 

ip,q) = ip',q')/{p" + q"), (19) 

the roles of A and C being interchanged after the transformation. 

\f A ^ = C, then a nonzero C is restored by a translation in the {p,q) plane. If 
A = 0, then the metric of the (p, g)-surface is brought back to the standard Szekeres form 
with y4 7^ 7^ C by a Haantjes transformation (a conformal symmetry transformation of a flat 
space, see [jl6|] for a description) in the {p,q) surface, which also restores the appropriate form 

of Qrr- 

The surface area of a (t = const, r =const) surface is finite only in the e = +1 case, 
where it equals A-nR^. In the other two cases, it is infinite. 

The 6 arbitrary functions /, M, a, P, Q and S represent 5 physical freedoms to control 
the inhomogeneity, plus a coordinate freedom to rescale r. 
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The density and Kretschmann scalar are functions of all four coordinates^ 

2(M' - 3ME'/E) 



4_2 8_ 2 

3P -gPP + Sp^ 



(21) 



where 

S^V^^ (22) 

is the mean density within "radius" r. For all p and p we have /C > 0, but assumptions of 
positive mass and density require p > and p > 0. Clearly there are density and curvature 
singularities at i? = — the bang and/or crunch — and at R' = RE' /E, M' ^ 3ME'/E — 
shell crossings. Additionally, p but not /C passes through zero where E' / E exceeds M'/3M. 

The matter flow = S", with projection tensor hafs = Qap + UaUj3, has the following 
properties, which are almost trivial to calculate with GRTensor 



" (R'-RE'/E) ' ^ ^ 

a" = M^V/3U" = 0, (24) 



= ^7°nV(,«,)+W(,a,))-f /i^ = ||^-||;^diag(0,2,-l,-l), (25) 

U;„/3 = V[/3Ua] + M[/3a„] = 0, (26) 

M(R' - RM'/3M) , , s , ^ 

Haf3 = \e^^^,C^^''psu<u' = 0. (28) 

Note that the relation between the active gravitational mass M and the "sum-of-rest- 
masses" M. is the same in the e = +1 Szekeres model as in the LT model: 



M' = M'/^l + /. (29) 

The sum of the rest masses contained inside the sphere of coordinate radius r at the time t is 
defined by: 

M = j p^dsx, (30) 

where (73 is the determinant of the metric of the t = const hypersurface, and the integral is 
taken with respect to the variables p and q from —00 to +00, and with respect to r from tq 
at the origin to the current value r. We have 

^ There are only two independent curvature invariants in the Szekeres metric, for which a good choice 
would be: R = Snp and C'^'^^Cap-yS = A8KP{(R' - RM' /iMf} / {R^{R' ~ RE' / Ef} = (4/3)(87rp)(87rp- 
87rp)^ — a pure Ricci invariant and a pure Weyl invariant. Though less tidy, the 87rp and K used above 
will also suffice. We thank a referee for pointing this out. 
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E 



VT+f 



{RlE)\RlEy. 



(31) 



Consequently 



1 r+co 

M = — dq 



dp / dx 



{t,p,q,x). 



(32) 



The term containing E' is integrated by parts with respect to x in order to move the prime 
(which, in the integrand, means ^) away from E to functions that do not depend on p and q. 
The result is 



M = - 



M 



M 



r] — 



(ro) 



— / dg / dpE- 

477 J-oo J-oo 



r+oo r+oo 

H / dq dp I dxE'^ 

StT J-oo J~oo 



ro 



3Mf' 



M' 



2(1 + ^/r+T, 



We note that 



(33) 



+ 00 



dg / dpE^"^ = Att 

J ~oo 

(this is the surface area of a unit sphere), and so 



(34) 



M 



M 



M 



VT+7 



vr+7 



(ro) 



+ 



3Mf' 



M' 



2(i + /)3/2 



dx. 



(35) 



From here, we obtain the same relation that holds in the L-T model, ( pOf) . 

Note that this result holds only in the e = +1 Szekeres model (the quasi-spherical one). 



With e = or e 
be defined. 



1, the total surface area of the (p, g)-surface is infinite, and so Ai cannot 



2.1 Special Cases and Limits 

The LemaTtre-Tolman (LT) model is the spherically symmetric special case e = +1, E' = 0. 
The vacuum case is (M' — 3ME' /E) = which gives M oc E^, and this requires 

M' = = S' = P' = Q' = E' (36) 



and any region over which this holds is the Schwarzschild metric in LT coordinates 
mass M. (See Jl^ for the full transformation in the general case.) 



with 



In the null limit, / ^ oo, in which the 'dust' particles move at light speed [|TT[ the 
metric becomes a pure radiation Robinson-Trautman metric of Petrov type D, as given in Exact 
Solutions [|T^, equation (24.60) with (24.62)[p'he Kinnersley rocket [0] is the e = +1 case of 



this null limit, which is actually more general than the axially symmetric form given in [|TI 



^Ref [|j corrected ref (T^'s mistaken claim that the null limit of Szekeres was a new metric. 
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The KS-type Szkeres metric was shown in [|TT] to be a special case of the above LT-type 
metric, under a suitable limit. 



2.2 Basic Physical Restrictions 

1. For a metric of Lorentzian signature ( — h ++), we require 

e + />0 (37) 

with equality only occuring where {R' — RE' / /{e + f) > 0. Clearly, pseudo-spherical 
foliations, e = —1, require / > 1, and so are only possible for hyperbolic spatial sections, 
/ > 0. Similarly, planar foliations, e = 0, are only possible for parabolic or hyperbolic 
spatial sections, / > 0, whereas spherical foliations are possible for all / > —1. 

2. We obviously choose the areal radius R to be positive, 

R>0 (38) 

(i? = is either an origin, or the bang or crunch. In no case is a continuation to negative 
R possible.) 

3. The mass M(r) must be positive, so that any vacuum exterior has positive Schwarzschild 
mass, 

M > 0. (39) 

4. We require the metric to be non degenerate &l non singular, except at the bang or crunch. 
Since (dp"^ + dq^)/E'^ maps to the unit sphere, plane or pseudo-sphere, \S{r)\ 7^ is 
needed for a sensible mapping, and so S > is a reasonable choice. In the cases e = 
or —1, E necessarily goes to zero at certain (p, g) values where the mapping is badly 
behaved. For a well behaved r coordinate, we do need to specify 

(R'-RE'/Ef 

i.e. (e + /) > except where {R! - RE'/Ef = 0. (41) 

In Lemaitre-Tolman models ^ {E' = 0, e = 1), the equality (1 + /) = = {R')^ 
can occur in closed models where the areal radius on a spatial section is at a maximum, 
or in wormhole models where the areal radius is minimum, R'{t,rm) = 0, V t. These 
can only occur at constant r and must hold for all (p, q) values. We will consider maxima 
and minima again later. 

5. The density must be positive, and the Kretschmann scalar must be finite, which adds 

either M' - 3ME'/E > and R' - RE'/E > (42) 

or M' - 3ME'/E < and R' - RE'/E < 0. (43) 

If (R' — RE'/E) passes through anywhere other than a regular extremum, we have 
a shell crossing, where an inner shell of matter passes through an outer shell, and the 
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density diverges and goes negative. This phenomenon is probably due to the spacetime 
coordinates being attached to the shells of matter, and is not physically realistic. Never- 
theless, we would like to avoid models in which such unphysical behaviour occurs, so it 
is useful to find restrictions on the arbitrary functions that prevent it. 

6. The various arbitrary functions should have sufficient continuity — and piecewise 
— except possibly at a spherical origin. 



3 The Significance of E 

3.1 Properties of E{r,p,q). 

Note that the Szekeres metric is covariant with the transformations r = g{f), where g is an 
arbitrary function. Hence, if i?' < in the neighbourhood of some value r = tq, we can 
take g = 1/f and obtain dR/df > 0. Therefore, R' > can always be assumed to hold in 
some neighbourhood of any r = tq. However, if R' changes sign somewhere, then this is a 
coordinate-independent property. 



As seen from eq. (|13D, with e = +1, E must be always nonzero. Since the sign of E is 
not defined by the metric, we can assume that E > 0. 

Can E' change sign? 

E' = {- [ip - Pf + {q- Qf] /S' + e}-\[{p- P)P' + {q- Q)Q']. (44) 



The discriminant of this with respect to (p — P) is 

1 



^2 



-^{q-Qf-'^^{q-Q)Q' + P" + eS'' 



The discriminant of Ap with respect to (g — Q) is 



(45) 



(46) 



Since, with e = +1, this is never negative, the equation £" = will always have at least one 
solution (exceptional situations), and in general two. The two exceptional situations are when 
Aq = 0. They are: 



1. S' = 0. Then E' = has a family of solutions anyway, but the solutions define a straight 



line in the (p, g)-plane. This will be dealt with below (see after eq. (|56| )) 



2. S' = P' = Q' = 0. Then i?' = at this particular value of r, and we see from eq. ([ 
that p will be spherically symmetric there. (In this case, the positions of the great circle 
from eq. ( [55D and of the poles from eq. (0) are undetermined). 



When Aq > 0, Ap will change sign at the following two values of q: 

S 



qi,2 = Q + gj {-Q' ± ^JP" + Q" + ^s") 



(47) 
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E' > E' > I E' > Delta p < 



q2 



q1 



E'>0 




E' > E' > 



E' > Delta p > 



E' > Delta p < 



p1 p2 



P ^ 



Figure 2: When S' > 0, E' < inside the circle defined by (^) and E' > outside. 
Note that the (p, (?) plane is a Riemann projection of a sphere, and on the sphere "inside" 
and "outside" are topologically equivalent. 



For every q such that qi < q < q2 there will be two values of p (and one value of p when q = qi 
or q = q2) such that E' = 0. Those values of p are 



q-Q Q'\' P'^ + Q'' 



S ^ S'J 



S'2 



(48) 



The regions where E' is positive and negative depend on the sign of S' . If S' > 0, then E' > 
for p < pi and for p > p2, if S' < 0, then E' > for pi < p < p2. E' = for p = pi and 
p = p2, but note that pi and p2 are members of a continuous family labelled by q. All the 
values of p and q from (|47| ) - ( |48D lie on the circle 



P-[P-P'§ 



+ 



Q-Q — 



P" + Q' 

5'2 



+ e 



The center of this circle is in the point 
and the radius of this circle is 



S' 



Ltp' 



E'=0 



'P'2 + Q 
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(49) 



(50) 



(51) 



The situation on the (p, g)-plane when 5" > is shown in Fig. |]. 



3.2 Properties of E{r, 9, 0). 

We consider the variation of E(r,p,q) around the spheres of constant t and r. 
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Setting e = +1 and applying the transformation ( |I6D to (0) and to its derivative gives 

E = r^^, (52) 
1 — cos 6 



, S' cos9 + sm9{P' cos(f) + Q' sini 

E — n 

1 — cos 



(53) 



E" 



S" cos 6 + sin e{P" cos + Q" sin <p) 
(1-cos^) 

/S'cos^ + sin^(P'cos0 + Q'sin< 
^ 1^^; 1^ (1-cos^) 

, {S'f + {P'f + {Q'f 



S 



(54) 



The locus i^^' = is 



S' cos 9 + P' sin 6* cos (p + Q' sin 6* sin = 0. (55) 

Writing z = cos6, y = sin ^ cose/), x = sin ^ sin 0, clearly puts {x,y,z) on a unit sphere 
through (0, 0, 0), and ( |55| ) becomes S'z + P'x + Q'y = which is the equation of an arbitrary 
plane through (0,0,0). Such planes all intersect the unit sphere along great circles, therefore 
i?' = is a great circle, with locus 

—5" 

tan0 = . (56) 

P' cos (p + Q'sm(p 



The plane has unit normal (P', Q', S')/^{P'y + (Q'^ + (S'^ . 

Now it is easy to understand the meaning of the special case S' = mentioned after eq. 
(|45|). As seen from (0), with S" = we have 9 = 0, which means that the great circle defined 
by P' = passes through the pole of Riemann projection. In this case, the image of the circle 
E' = on the {p,q) plane is a straight line passing through {p,q) = {P,Q), as indeed follows 
from (|44|). The sign of E' is different on each side of the straight line. Compare also with Figs. 
I and §. 



From ( |53|) and (|52| ) we find 

E' S' cos e + sme{P' COS + Q' sin 



(57) 

E' 



E S 

thus 

constant =^ S' z + P'x + Q'y = S x constant (58) 



E 

which is a plane parallel to the E' = plane, implying that all loci E' / E = constant are small 
circles parallel to the E' = great circle. This will be seen to apply to shell crossings and 
apparent horizons. 

The location of the extrema of E' / E are found as follows 

djE'/E) sm9{P'sm<p-Q'cos<P) 

= s = ° ^ (^^^ 
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Q' P' 
tan^e = 7^ ^ cos0e = ei , ei = ±1, (60) 

= ^ (61) 



d{E'/E) S' sin 9 - P' cos cos - Q' cos 9 sin , 



^ ^ p-cos0e + g^sin0e v/(^o^ + (Q'y ^ .... 

tan^e = = ei^ (62) 

S' 

COs9e = 62 , £2 = ±1. (63) 

^{sr + {P'Y + {Q'Y 



The extreme value is then 



-e] . =-'^- s • 



Since (sin^ecos0e,sin^esin0e,cos^e) = t2{P\Q',S')l^{P'Y + {Q'Y + {S'Y , eq. (^ 
shows that the extreme values of E' /E are poles to the great circles of E' = 0. The latter can 
now be written in parametric form as 

cos 6^ = — cos ^/^ sin 6'e , (65) 

cos 6'e tan 0e + tan '?/' 

tan0= ■ (66) 

cos o'e — tan 0e tan ip 

Clearly E' / E has a dipole variation around each constant r sphere, changing sign when 
we go over to the anitipodal point: {9, (p) ^ (n — 9,(f) + n). Writing 

we see that RE'/E is the correction to the radial separation R' of neighbouring constant 
r shells, due to their not being concentric. In partuicular RS' / S is the forward {9 = 0) 
displacement, and RP'/S &l RQ' / S are the two sideways displacements (9 = n/2, (p = 0) &l 
[9 = n/2, = 7r/2). The shortest 'radial' distance is where E' / E is maximum. 



It will be shown in section |]2| that, where R' > 0, E'/E < M'/{3M) and E'/E < R'/R 
are required to avoid shell crossings, and also in eq ( |130|) that R'/R > M'/3M. These 
inequalities, together with M' > 0, imply that the density given by (^), as a function of 

X = E'/E: 



2M' 1 - 3Mx/M' , , 

l-Rx/R' 



has a negative derivative by x: 



_R/R! - 3M/M' 2M' 

{1-Rx/R'y 'Wr^^' ^ ^ 
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and so the density is minimum where E' / E is maximum. 



The density, eq. can be decomposed into a spherical part and a dipole-like part, as 

noted by Szekeres [|5] and de Souza |P (see also p. 30 in Ref.|p.7[). Rewriting de Souza's 
result into our notation of eqs. (|l|) and (^, we obtain: 



p = Ps + Ap, (70) 



where ps is the spherical part: 



_ 2M'{A + C) - m{A' + C) 

~ R^[R'{A + C)-R{A' + C')y ^ ' 

and Ap is the dipole-like part: 

_ A' + C' -{A + C)E'/E 6MR' - 2M'R 

^~ R' - RE'/E R^[R'{A + C) - R{A' + C')]' ^ ' 

The dipole-like part changes sign on the surface where E' / E = {A' + C')/{A + C), but lacks 
the antisymmetry property: Ap{—E'/E) ^ —Ap{E'/E). It can be verified (see Appendix 0) 
that the Ap = hypersurface does intersect every (t = const, r = const) sphere along a circle, 
unless P' = Q' = S' = {= A' = C), in which case the dipole component of density is simply 
zero. The surface Ap = in a t = const space is comoving, i.e. its definition does not depend 
on t. Also, its intersection with any sphere of constant r, E' / E = {A' + C) /{A + C) = const, 
is a circle parallel to the great circle E' = 0, as noted after eq. (^8]). It will coincide with the 
E' = circle in those points where A' + C = (if they exist). The dipole-like component 
will be antisymmetric with respect to E' / E only at such values of r, where {A + C)R = 
= {A' + C')R', but such values may exist only at the center, R = 0, because A + C = 
contradicts eq. (||). 

In the maximum (e2 = —1) and minimum (e2 = 1) directions, 

F" = 

max/min 



S{S"S' + P"P' + Q"Q') - [{S'f + {P'f + (g')'](25' + e2J{S'Y + {P'f + {Q'f 



(73) 



S{S' - e2^[S'Y + {P'? + {Q'? ) 

while around the E' = Q circle 

^„ {S'f + {P'f + {Q'f 5"(P'cos0 + Q'sin0)-y(P"cos0 + Q"sin0) 

E = / • (74j 

S J {P' COS (j) + Q' sin (Py + {S'y - (P' cos + Q' sin 0) 



4 Regular Origins 

When e = +1, R = occurs at an origin of spherical coordinates, e.g. R(t,0) = 0, V t, 
where the 2-spheres have no size. Similarly, R{t,0) = = R{t,0), etc V t. There will be a 
second origin, at r = rp say, in any closed, regular, / < model. Thus, by (|Tl|) and (^ and 
their combinations with (^) &l (^), for each constant r] 

M P 
lim — = , lim / = , lim — = 0. (75) 
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The type of time evolution at the origin must be the same as its neighbourhood, i.e., along 
a constant t slice away from the bang or crunch, by (0) and (0), 



< lim 



\f\y\t 



Clearly, we need M — 0, / — and 



< lim 



M 

f|3/2 



< oo. 



Using I'Hopital's rule, this gives 



lim 



M 
3Mf' 



-0 2M'f 



< oo. 



1. 



(76) 



(77) 



(78) 



The density and Kretschmann scalar must be well behaved. We don't consider a vacuum 
region of finite size at the origin, as that is just Minkowski space, M = 0, and we don't consider 
the obscure case of a single vacuum point at the origin. Because p (Si p in (p]) evolve differently, 
we also need 

, 6M 2M' 3R'M ^ 

< lim = lim ^^tttt: < oo =^ lim ^rrrr = 1 (79) 



^0 R3 



lim „ 



r^O RM' 



and 



, 2(M' - 3ME'/E) , 2M' (1 - 3ME'/M'E) 

D < lim — — - = lim < oo 

R2(R' - RE'/E) r" Rm' (1 - RE'/R'E) ' 



(80) 



but in fact the latter is ensured by the former, and the anisotropic effect of E vanishes at the 
origin. However, since E' / E is restricted by the conditions for no shell crossings, it would be 
odd if lim^^o ^-£^7^'-^ oi" iii^r-^o RE' / R' E were divergent. Since 

R' ( 3Mf'\ 

Mr 



M' 



+ 



M2 



-ff 



sin ri{ri — sin r^) (1 — cos r]) 
cosr])^ 



(81) 



eqs. (^) and (|7q ) above make the first term zero and the second non-zero at an origin for all 

< 7] < 2n, so we only need 

Ma' 

hm < oo. (82) 

r^O M' 

Lastly, the metric must be well behaved, so E should have no unusual behaviour, such as 
S = 0, that would compromise a valid mapping of (dp"^ + dq^)/E'^ to the unit sphere. Also, 
to ensure the rate of change of proper radius with respect to areal radius is that of an origin, 
grr/{R'Y should be finite 



^ {R'-RE'/Ef 



< lim 1 



and 



lim 



3ME'^ 
M'E 

ME' 
M'E 



< oo 



oo < lim 



ME' 



M'E 



< oo 



(83) 
(84) 
(85) 
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where the last of ([79|) has been used. This should hold for all (p, g), i.e. all 
gives 



Thus (|7]) 



oo < lim 



MS' 



M'S 



< oo , 



oo < lim 



MP' 



M'S 



< oo 



oo < lim 



MQ' 



M'S 



< oo, (86) 



all three limits being different from 1/3. 

All of the above suggests that, near an origin, 

M ~ , / ~ , 5 ~ i?" , P ~ i?" , Q ~ i?" , n>0. 



(87) 



The condition E' / E < M' /2>M that will be obtained in the next section implies that near an 
origin 

n < 1 (88) 



5 Shell Crossings 



5.1 Occurrence and Position of Shell Crossings in a Surface of Con- 
stant t and r. 



A shell crossing, if it exists, is the locus of zeros of the function R' — RE' /E. Suppose that 
R' - RE'/E = holds for all r at some t = to. This leads to S' = P' = Q' = R' = 0. Since 
P, Q and S depend only on r, this means they are constant throughout the spacetime. As 
seen from (P and (^Of), the Szekeres metric reduces then to the LT metric, and so this case 
need not be considered. 

Suppose that R' — RE' /E = holds for all t at some r = tq. This is an algebraic equation 
in p and q whose coefficients depend on t and r. Taking the coefficients of different powers of 
p and q we find P' = Q' = S' = R' = 0, but this time these functions vanish only at r = tq, 
while R'{t,ro) will vanish for all t. This will either be a singularity (when M'(ro) 7^ 0) or a 
neck (when M'(ro) = 0), familiar from the studies of the LT model, see Refs. [jlO|] and [|13| . 
Hence, R' — RE'/E ^ except at a shell crossing or at special locations. 

Now R' > and R' - RE'/E < cannot hold for all p and q. This would lead to 
E' > ER'/R > 0, and we know that E' cannot be positive at all p and q. Hence, with 
R' > 0, there must be a region in which R' — RE'/E > 0. By a similar argument, R' < and 
R' — RE'/E > cannot hold for all p and q, so with R' < 0, there must be a region in which 

R' - RE'/E < 0. 

Assuming R' > 0, can R' — RE'/E be positive for all p and g? Writing 



ER'/R - E' 



1 

2S 



'S^ R^ 



1 



-eS 



S 



IT 
R 



S' 



P)P'+(q-Q)Q'], 
the discriminants of this with respect to {p — P) and {q — Q) are 

p/2 I 



^ . o rt 



l + f ) (q ' Qf + {q - Q)Q' 



/S' 






-5)1 



(89) 



(90) 
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_ 1 /s; IT 



52 



(91) 



Thus ER'/R — E' will have the same sign for all p and q when < (because then also 
Ap < for all q). Hence, ER'/R — E' has the same sign for all p and q (i.e. there are no shell 
crossings) if and only if 

> ^ ^2 ■■= (^)- (92) 

Note that when e = 0, this can fail only at those points where R' = 0. 

If R!"^ / R^ = $2, then A^ = 0, and so Ap = at just one value of g = qss- At this value 
of q, ER'/R — £" = at one value of p = pss- In this case, the shell crossing is a single point 
in the constant (t, r)-surface, i.e. a curve in a space of constant t and a 2-surface in spacetime. 

If R''^/R^ < $2, then the locus of ER'/R-E' = is in general a circle (a straight line in 
the special case S'/S = —R'/R) in the (p, g) plane. The straight line is just a projection onto 
the (p, g) plane of a circle on the sphere of constant t and r, and so is not really any special 
case. 

When Ag > {R'^ /R^ < $^), the two limiting values of q at which Ap changes sign are 



-Q' ± + g'2 _^ e (5/2 _ sm'^R^) 
^^'^ = S'/S + R'/R ' ^^^^ 

and then for every q such that qi < q < q2, there are two values of p (only one \f q = qi or 
q = ^2) such that ER'/R, - E' = 0. These are 

-P' ± J- [(f + f ) (g - Q) - Q'] ' + P'2 + g/2 + , (5.2 _ sm'^/R^) 

= S'/S + R'/R • 

The values of p and q from (153) and (P?]) lie on the circle with the center at 



( P' Q' \ 

{psc, qsc) =[P- s'/s + R'/R^ ^ - S'/S + R'/R ) ' (^^^ 



and with the radius Lsc given by 



2 + Q'^ + e {S'^ - S^R'yP^) 



- {S'/S + R'/Rf • 



This is in general a different circle than the one defined by E' = 0. As seen from (|89| ), the 
shell crossing set intersects with the surface of constant t and r along the line E' / E = R'/R = 
const. As noted after eq. (pB]), this line is a circle that lies in a plane parallel to the E' = 
great circle. It follows immediately that the E' = and the SC circles cannot intersect unless 
they coincide. 



5.2 Conditions for No Shell Crossings 

Szekeres p3] obtained a number of regularity conditions for the e = +1 metric, namely: (1) 



On any constant time slice, R(t = const, r) is monotonic in r, which allows a transformation 
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to make R = r and _R' = 1 on that slice. (2) At an origin, A, Bi, B2 &i C should be C^, 
/ = 0, and M but we are not sure why he required A' = Q = B[ = B'^ = C there. (3) 

To keep the density non-singular, < {{S'f + {P'f + {Q'f)/S'^ < min{{R'lRf, {M' /?,Mf), 
which is a no shell crossing condition. We shall improve on the latter below. 

For positive density, (|0]) shows that (M' - 3ME'/E) L {R' - RE'/E) must have the 
same sign. We now consider the case where both are positive. Where (M' — 3ME'/E) < 
and (R' — RE'/E) < we reverse the inequalities in all the following. 

In the case of both (M' - 3ME'/E) ^ {R' - RE'/E) being zero, this can hold for a 
particular (p, g) value if M' /"iM = R' / R, but the latter cannot hold for all time. This case 
can only hold for all {p,q) if M' = 0, E' = 0, R' = 0, which requires all of M', /', a', S', P' , 
Q' to be zero at some r value. 

We consider the inequality (Af — 3ME' /E) > and we argue that it must hold even for 
the extreme value of E' / E, (|54D, for which we obtain 

= ^ V r. (97) 

c 

It is obvious that this is sufficient, and also that 

M' > V r. (98) 

We will now consider {R! — RE'/E) > for all 3 types of evolution. 

5.2.1 Hyperbolic evolution, / > 

For hyperbolic models, we can write: 

R' M' . ^. /'/3, \ P'^a' 



M' E' 
> 



3M - E 



where 



At early times. 



(1 - 04) + ^ (-04 - 1) - L^05, (99) 

(100) 



R ^ ' f \2 

sinh?7(sinh?7 — r]) sinh?7 



(cosh?7 — 1)2 ' (cosh?7 — 1)2 



(101) 



R - y^ + O(r^^)-0, (102) 
05 ^ 4 + (103) 



we find ^5 dominates and 



'r ^ iW"*^^' 



(105) 
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so that {R' - RE'/E) > gives 

a' < V r. (106) 

Similarly, at late times, 

7] ^ oo , _R— >oo, 05— >0, 04— (107) 
we find 05 vanishes and 

- - -4, (108) 
R 2 f ^ ^ 

so that 

Following the above analysis of (Af — 3ME'/E) > we obtain 



Yf> s (110) 

which obviously implies 

/'>0 Vr. (Ill) 
Again, since we already have M' > 0, it is clear that this is sufficient, and that 

R' > 0. (112) 

5.2.2 Parabolic evolution, / = 

The easiest way to obtain the conditions for this case, / = 0, /' 7^ 0, is to put f] = rij \ff > 
in the hyperbolic case, and take the limit f ^ 0, r] ^ 0. All terms involving /'// cancel and 
we retain exactly the same conditions, viz ( p.06| ) &i ( |lllj ) (and of course (|97|) ). Naturally, ( |110|) 
ceases to impose any limit. 



5.2.3 Elliptic evolution, / < 

For elliptic models, we can write: 



where 



At early times. 



^ sin 77(77 — sin 77) ^ sin 77 

{l — cosr)Y ' (1 — cos77)^ 



(113) 
(114) 



7/ ^ 0, (115) 

M 

R - (Z7)y + ^(^')^0' (116) 

02 ^ ^ + 0(77)^+00, (117) 

01 - ^+^(^')^| (118) 



we find 02 dominates and 

R ^ 

so that R{R'/R - E'/E) > gives 



(119) 



a' <0 V r. (120) 

Similarly, at late times, 

V ^ 27r, (121) 

R - (^^^ + O((2.-,)^)-.0, (122) 

02 ^ -4/(27r - 77)3 + 0((27r - r/)) ^ -oo, (123) 

01 ^ -87r/(27r-r/)3 + 2/3 + 0((27r-r7)) ^ -oo, (124) 



we find 



^ M' ( Svr \ f / 127r \ (-/)3/^a- / 4 \ 

i? M^27r-r/)3j / ^27r - r/)^ j ^ M [{27r - t]^ J ' ^ ' 

so that R^I^{R'/R - E'/E) > now gives 

27rM /M' 3/'\ , , , , 

' — - ^ +a'>0 Vr, (126) 



{-fr/'\M 2/ 

which is the condition that the crunch time must increase with r. Since we already have 
M' > 0, it may be easily verified that these conditions are sufficient to keep 

R' >0 (127) 

for all r]. 

We now show the above also ensure R(R' / R — E'/E) > always. Defining the crunch 
time 6(r) with 

27rM 

^ = «+(Z7)^ (128) 

we can re-write ( |113|) as 

R' M' b' [2 \ i-a') [2 ^ ^ \ 

+ 71 T o - + 71 T o - + 2vr02 . (129) 



R 3M ih-a) V3 V (6 - a) V3 

The derivative of (2/3 — 0i) is [2,7] — 3 sinr/ + cos ?7)/(l — cos r/)^, and the third derivative of 
the numerator of the latter is rj sinr^. It follows that (2/3 — 0i) > and declines monotonically 
from +00 to as goes from 27r to 0. Since (2/3 — 0i + 27r02) is the mirror image \n rj = n 
of (2/3 — 0i), we have that 

R' M' ^^^^^ 

R > 3M' 

so that (|97| ) guarantees that for each given r, the maximum of i?'/_E' as {p,q) are varied is no 
more than the minimum of R' / R as t] varies. 
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Note that although (|126|) implies 

f M' 

a condition such as (|110| ) is not needed in this case. As an indication of the approximate 
magnitude of R' /R\ram, at the moment of maximum expansion along any given worldline, 

— = — , (132) 

so it would be possible to have E' / E\„ia.x close to i?'/-R|min around the time of maximum 
expansion. 



6 Regular Maxima &. Minima 

Certain topologies necessarily have extrema in R. For example, closed spatial sections have a 
maximum areal radius, and wormholes have a minimum areal radius, i.e. R'{t,rm) = 0, V t. 

Suppose (e + /) = at some r = r^- By we must have 

firm) = (133) 

(unless /' is discontinuous there, which we won't consider). We need {R' — RE' / E) = to 
keep Qrr finite, and hence (Af — 3ME'/E) = to keep p finite, both holding V {t,p,q) at 
that Tm- More specifically, along any given spatial slice away from the bang or crunch, we want 

{R' - RE'/E) 



{M' - 3ME'/E) 



> L , 0< L <oo, (134) 
N , 0<N<oo. (135) 



{R' - RE'/E) 

As noted above, we require 

M' = f = a' = S' = P' = Q' = (136) 

to ensure 

R' = 0. (137) 
The limits ( p4l ) and ( |T35| ) must hold good for all t and for all {p,q), so using (|13D, (0) &i 



([IT3| ) with R> 0, M > 0, S > shows that 

M' f a' R' 



S' P' Q' E' 



(138) 
(139) 



must all have finite limits, that do not have to be zero. Using I'Hopital's rule, each of the 
above limits can be expressed in the form 

M' 2M"^F+7 2M" 

= M 7^ = — T — = -Li- 
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Thus, for / = —1, the above condidtions for no shell crossings in elliptic regions should be 
re-expressed in terms of these limits. 

It is worth pointing out that E' = at f = —I does not imply the shells near an extremum 
in R are concentric. It is the above limits that determine whether there is non-concentricity at 

Conversely, imposing R' = forces all of (|136|), if we are to avoid shell crossings. To 



obtain / = — e, we must impose one further requirement for a regular extremum — that no 
surface layers should occur at r = r^. Using the results for the normal n^j and the extrinsic 
curvature Kij shown in the next section, and choosing the junction surface to be at constant 
coordinate radius, r = Z = r^, the non-zero components are: 

n. ^ JJ^-^, (141) 

K„ - (143) 

Now at an extremum in R, the factor {R' — RE' / E) goes from positive to negative (because 
where {R' — RE' /E) < 0, the no-shell-crossing conditions require R' < 0), which means that 
flips direction. For a boundary with no surface layer we must have pointing the same 
way on both sides, towards increasing r say, and zero jump in the extrinsic curvature. So, if we 
cut the model at the maximum or minimum r,^, and match the two halves back together, we 
need 

K±i-n,) = K-^i+n,) K± = -K^, (144) 

which is only possible if 

/ = -e. (145) 

If however, is only a shoulder — i.e. R'{rm) = 0, but R' has the same sign on either 
side, then the normal direction does not change sign, so there is no surface layer even if / 7^ — e. 
However grr = L"^ goes to zero, so it is likely that a change of coordinates could make \R'\ > 0. 



6.1 Summary: Conditions for No Shell Crossings or Surface Layers 

The conditions found here are exactly those on M, f &l a for LT models (see which 
generalises those of [|| for a = LT models), with extra conditions involving S, P, Q also. 
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e 


R' 


/ 


M', f, 


a' 




S', P', Q' 


+1 


> 


all 


M' >0 


v/(5')2+(P')" + (Q')2 ^ M' 

o — Oik/ 


+1 


> 


> 


f'>0 
a ^ u 

but not all 3 equalities at once 


S - 2/ 

(no condition where / = 0) 


+ 1 


> 


< 


2nM ( M' 3f'\ 1 „/ ^ n 
(-/)3/2 \M 2/ y " — " 

a' < 

but not all 3 equalities at once 




+ 1 


= 

R" > 
neck 


-1 


M' = 
(/ = - 

2nM 


, /' = 0, a' = 

1 for no surface layer) 


s' = 0, p' = 0, g' = 

v/(5")2+(P")2+(Q")2 ^ M" 


a" < 


S — 3M 


+1 


= 

R" < 
belly 


-1 


M' = 
(/ = - 

2itM 


, f = 0, a' = 

1 for no surface layer) 

(^-f)+a"<0 


s' = 0, p' = 0, g' = 


a" > 


5 — 3M 


+ 1 


< 


all 


M' < 


+ + ^ A/' 


+ 1 


< 


> 


/'<o 

a' > 

but not all 3 equalities at once 


V(S')2+(P')2+(Q')2 ^ ^' 
5 - 2/ 

(no condition where / = 0) 


+1 


< 


< 


a' > 

but not all 3 equalities at once 





7 Impossibility of a Handle Topology 

Since the function E has the effect of making the distance between adjacent constant r shells 
depend on angle, this allows us to create a wormhole that is bent, so that the two asymptotic 
world sheets on either side can be thought of as intersecting in the embedding. 

This leads to the question of whether those two world sheets can be smoothly joined across 
a junction surface. In fact the possibility of matching the two world sheets together can be 
considered independently of whether there is a natural embedding that would allow them to 
intersect at the appropriate angles. 

Thus we investigate whether it is possible to create a Szekeres model with a handle topology 
in the following way. Take a wormhole topology — an e = +1 model with r = 0, /(O) = — 1 
at the wormhole &l f < Q nearby — and let it be mirror symmetric about r = 0. Choose a 
comoving open surface S on one side of the wormhole, and its mirror image on the other side, 
and match the two sheets together along S, as shown schematically in Fig. |. Because the 2 
sheets are mirror images, this is equivalent to matching S to its own mirror image. 
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Figure 3: Conceptual illustration of joining a Szekeres wormhole model to itself across a 
boundary surface E, shown as a heavy line. The boundary may be close to the wormhole, as 
shown, or out in the asymptotic regions. There is no significance to the change from solid to 
dotted circles, other than picture clarity. The handle topology is shown as an embedding of 
a constant time section, with one angular coordinate suppressed, although a valid matching 
across S does not require the result to have a natural embedding. However, it is shown that 
the matching fails because it is not preserved by the model evolution. 



To implement this we choose a comoving surface: 

rj: = Z{p,q) (146) 

and surface coordinates: 

e = it,p,q). (147) 

The two fundamental forms and the normal are calculated in appendix 

Obviously 1st fundamental forms match by construction, and normal vectors are equal and 
opposite: 

n+ = -n-. (148) 

In fact the 2 surfaces S+ &l S_ are identical except for the sign of n^. Thus K^j = —Kij, so 
the only way to match the 2nd fundamental forms, = K~j, is to make them zero: 

= 0. (149) 

The Kpt = —Kpt &i Kgt = —Kqt equations give 

R' = =^ rs is at a shell crossing unless Z = const (150) 
or i? = Static: Not possible (151) 

or Zp = , Zg = rY. = Z = constant. (152) 
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If the matching surface is at constant r, then only a closed torus topology is possible. So the 
answer is: no, a handle topology is not possible. 

Suppose R = possible, then it should be possible to solve 

Kpp = , Kp, = , Kgg = (153) 

for Z{p,q), by specifying suitable functions E{r,p,q) &l R(r) on an initial time slice. In other 
words, you can probably match on a constant time slice, but the matching is not preserved by 
the model evolution. 



8 Szekeres Wormholes? 



it has been shown in []I0| that LT models can describe the Schwarzschild-Kruskal-Szekeres 
manifold, as well as models that have the same topology but non-zero density. It has also 
been shown that the matter flows from past to future singularity, with possibly some matter 
escaping to ^7+ or some being captured from J'^ . The effect of the introduction of matter on 
the causal structure is to split the Kruskal event horizons and reduce communication through 
the wormhole. The locus R = 2M is an apparent horizon, but not an event horizon, and light 
rays fall irrevocably through the AH towards the singularity wherever M' > 0. Only if the 
density is (locally) zero is i? = 2M (locally) null. Only if the density is everywhere zero is 
R = 2M the event horizon. (See also for a study of light rays and AHs in a collapsing LT 
model with a = 0.) 

Since LT models are a subset of Szekeres models, it is of interest to look at the properties 
of the Szekeres generalisation, and determine how the loss of spherical symmetry in the Szekeres 
model affects the LT result. 

In particular, given the anisotropy of the metric and the fact that the proper separation 
of constant r shells varies with p &i q, \s it possible for null or timelike paths to pass through 
a neck or wormhole, by choosing a path along which distances have been made shorter by the 
particular form of El In other words, can one construct a Szekeres wormhole that is traversible? 

For a wormhole, we require an elliptic region, in order to create a "neck" - a regular 
minimum in R(t = const, r), 

-l</<0 , e = +l, (154) 
but the asymptotic regions may be described by elliptic, parabolic, or hyperbolic regions. 



8.1 The Fastest Way Out 

The general null condition gives 



o = k-k^g^, = {-i){kr + ^^j^f^ikr+^,m'+{k'f) (155) 



e + / \dt J \ \dt J \dt 
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It is obvious that at each event dr/dt is maximised by choosing kP = Q = k'^. Since R is 
independent of {p,q), this also gives the direction of maximum dR/dt\nuU at any event. We 
will call this "radial" motion, and radial null paths "rays". Thus, the DE 

, dt 
f' = — 

" dr 

in principle solves to give 

t = tn{r) (158) 

along the " ray" . We don't expect this to be geodesic, but we regard it as the limit of a sequence 
of accelerating timelike paths, and thus the boundary to possible motion through a wormhole. 
The acceleration of this path may be calculated from a" = k^W pk"' , as given in appendix |C|. 



8.2 Apparent Horizons 

The areal radius along a 'ray' is 



Rn = R{tn{r),r), (159) 
{RJ = Rt'^ + R! (160) 

■ ^ V-^Ui?' (161) 



v/r+7 V E 



¥ + / RE' 



R'--—\+R!, (i = ±l. (162) 



' v/r+7 V E 

These rays are momentarily stationary when 

{Rn)' = 0. (163) 

Now light rays initially along constant p and q will not remain so, owing to the anisotropy of the 
model. However, since these "radial" directions are at each point the fastest possible escape 
route, we define this locus to be the apparent horizon (AH). 

(Szekeres ||2^ defined a trapped surface as the locus where null geodesies that are (mo- 



mentarily) 'radial' have zero divergence, k^-^ = 0, where k^-yk^ = 0, k^j,k^ = 0, A;^ = = k'^ . 
He obtained 

Given the anisotropy of the model, one doesn't expect this to be the same locus as our AH.) 

Assuming a normal spacetime point will have non-zero metric components, and taking R 
increasing with r on constant t slices, 

CRE'\ 
R! - -^J > 0, (165) 

we require 

(ij = -1, (166) 
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i.e. 



Either (future AH: AH+) j = 

i = 

Or (past AH: AH") j = 

i = 

Note that we want 'outgoing' to mean moving 
through the neck would change from incoming ■ 
there, j would also have to flip there. 



+ 1 (outgoing rays) 

—1 (in a collapsing phase) (167) 

— 1 (incoming rays) 

+1 (in an expanding phase). (168) 

away from the neck at r = 0. A ray passing 
to outgoing at r = 0, and, since R' flips sign 



8.2.1 The Apparent Horizon and its Location with Respect to E' = 0. 

Define 



D 



(169) 



Then 



{D>Q)^{R> 2M). (170) 
Since M/R and {2M / R + f) are positive, we see that D >1 leads to a contradiction, and so 

D<1. (171) 

However, \D\ can be greater than 1 because L> < — 1 is not prohibited. We have 

M \ 



{D < -1) 



R< 



1 + ,/1 + fJ- 

This will always occur when R is close to the Big Bang/Big Crunch. 
Using D, the equation of the AH is 

RE' + DR'E = 0, 

and in terms of p and q this equation is 



(172) 



S R 



S' ^R'' 



{p-Py + iq-Qr +2[{p-P)P'+{q-Q)Q']-S'{^ + D 



R 



The discriminant of this with respect to p is 



(173) 



= 0. 
(174) 



R 

The discriminant of this with respect to q is 



Ag = 64 



S R 



12 



12^ 



,R^ 



(175) 



(176) 
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Now, if Ag < everywhere, then Ap < for all q, in which case there is no p obeying 
(|174| ), i.e. the apparent horizon does not intersect this particular surface of constant 



If Ag = 0, then Ap < for all q except one value q = q^, at which Ap = 0. At this value 
of j9 = pq, (|174| ) has a solution, and so the intersection of the apparent horizon with this one 
constant (t, r) surface is a single point. Note that the situation when the apparent horizon 
touches the whole 3-dimensional t = const hypersurface at a certain value of t is exceptional, 
this requires, from ( |174|) , that P' = Q' = S' = R' = at this value of t. The first three 



functions being zero mean just spherical symmetry, but the fourth one defines a special location, 
as mentioned at the beginning of sec. [ST^. These equations hold in the Datt-Ruban [^, ^ 



solution. 

If Aq > 0, then Ap > for every q such that qi < q < q2, where 



-Q' ± + g/2 ^s^[^-D^^ 

5 

and then a solution of ( |174|) exists given by 



^^'^ - S'/S - DR'/R 



-P'± 



'^-D§){q-Q) + Q'] +P'^ + Q'2 + S^{^-D^^ 



P''' = S'/S - DR'/R • (^^^^ 

Except for the special case when S'/S = DR'/R, these values lie on a circle in the {p,q) plane, 
with the center at 

( P' Q' \ 

^^^^ S'/S-DR'/W ^ ' S'/S - DR'/R ) ' ^^^^^ 



and with the radius Lah given by 



P'^ + Q'^ + S^ (|^_D2H'^ 



- {S'/S-DR'/Rf 

The special case S'/S = DR'/R (when the locus of AH in the {p,q) plane is a straight line) 
is again an artefact of the Riemann projection because this straight line is an image of a circle 
on the sphere. 

In summary, the intersection of AH with the (p, g)-plane is 

• nonexistent when R''^ / R^ > (this is the same $ as for the shell crossing); 

• a single point when R'"^ / R^ = <I>^/D^; 

• a circle or a straight line when R'"^ / R^ < 



The condition R'"^ / R^ < <I>^/D^ is consistent with the condition for no shell crossings, eq. 
(|5^), when \D\ < 1. We already know that necessarily D <1, but D < — 1 is not excluded. 

With \D\ < 1, when the intersection of AH with {t = const, r = const) is a single point, 
a shell crossing is automatically excluded. 
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Figure 4: The circles CI and C2 on a sphere (seen here edge on) will project onto the 
plane (seen here as the horizontal line) as the circles PI and P2 that are outside each other. 
Only parts of PI and P2 arc shown here. Circle CI is the E' — set, circle C2 is the 
apparent horizon circle. 



Note that from ( |173| ) and from the assumptions R > 0, E > and i?' > we have 

(D > 0) ^ {E' < 0) 

(D < 0) ^ {E' > 0). (181) 

But D > and D < define regions independent of p and q. Hence, on that surface, on which 
D > 0, E' < on the whole of AH. Where D < 0, E' > on the whole of AH. This implies 
that the E' = circle and the AH cannot intersect unless they coincide. Indeed, these circles 
lie in parallel planes, by the same argument that was used at the end of sec. p] : the line on 
the (if:, r) = const surface defined by ( |173| ) has the property E' / E = —DR'/R = const, and 
so it must be a circle in a plane parallel to the E' = great circle. It follows that of the three 
circles {E' = 0, SC and AH), no two can intersect unless they coincide. 

When the E' = and AH circles are disjoint, they may either be one inside the other or 
each one outside the other. However, when projected back onto the sphere, these two situations 
turn out to be topologically equivalent: depending on the position of the point of projection, 
the same two circles may project onto the plane either as one circle inside the other or as two 
separate circles, see Figs. ^ and ^. 



8.2.2 Location of the AH Compared with R = 2M 

Along R = 2M 



(R^y = R\l + ij)-^, (182) 



so R = 2M is not the AH except where E' = 0. 
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Figure 5: The same circles as in Fig. ^ projected onto a plane from a different pole will 
project as one inside the other. The transition from the situation of Fig. ^ to that of Fig. 
H is continuous and occurs when the sphere is rotated, but the pole and the plane are not 
moved. Then one of the circles (CI when a clockwise rotation is applied to Fig. ^ will 
pass through the pole at one value = (/Jq of the rotation angle. Its image on the plane is 
acquiring a larger and larger radius with increasing until it becomes a straight line when 
Lp — LfQ. When Lp increases further, the straight line bends in the opposite direction so that 
it surrounds the second circle P2. 



Eq ( p.62| ) with ( |163|) and £j = — 1 can be written 

_ 2M(1 - RE'/R'Ef 

" l + f\l~{l-RE'/R'Ey] ^^^^^ 
- V^RE'/I^E. (184) 



1 + /(2\/ - 1/2) 



The effect of E(r,p,q) is to create a dipole in the geometry and density around each {t,r) 
shell, with E' = on an "equator", and extreme values 



E' 

'e 

at the poles. 



extreme 



(185) 



It is clear then that "radial" displacements between two nearby surfaces of constant r are 
shortest where E' / E is maximum, and light rays move outwards fastest (max dr/dt, min dt/dr, 
i.e. most rapid transfer between constant r shells at the same {p,q) value). It has also been 
shown that the density is minimum here. The longest "radial" displacements, slowest light ray 
motion and maximum density occur at the opposite pole. 

We will call the direction where E' / E is maximum, the 'fast' pole, and where E' / E is 
minimum, the 'slow' pole. 



Now the conditions for no shell crossings require 

E' 

'e 



M' 

3M ^ ' 

extreme 
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and for an elliptic region we have 

SO 

V^<1 , (l-\/)^>0 , -3 < (2\/-y2) < 1. (188) 

In those places on the AH, where = = E' , we see that the surface R = 2M intersects 
the AH at all times, but the AH is a kind of oval with half inside R = 2M and half outside. 

For / = -1 

^ = 1 (189) 

regardless of V. So it is clear that AH+ &l AH~ cross in a 2-sphere at the neck of the wormhole 
(/ = —1) at the moment of maximum expansion {R = 2M), as in LT. Note also that at 
the bang, wherever a' ^ 0, R ^ and R' ^ oo imply V ^ 0, and the anisotropy becomes 
negligible. Similarly for the crunch. 

But in general, for all > / > — 1 

^ > ^ > (190) 



1 + 31/1 - 2M 

and Rah/^M decreases monotonically as V goes from —1 to 1. Note that iKxtremel is likely 
to be less than 1, and also that the maximum &l minimum values of Rah do not have a simple 
relationship. 

We have that Rah/'^M < 1 where V > i.e. where E' > 0. In other words, taking 
a (t,r) shell that intersects the AH at the fast pole, the light rays move fastest between the 
shells exactly where the shell is just emerging from the AH. 

Some other features of the AH locus are discussed in appendix |D[ 



8.3 Causal Structure of A Szekeres Wormhole 

We shall next establish whether a radial null ray can pass through a Szekeres wormhole. We 
shall have to treat the neck separately from every other r value, because of the need to treat 
the / — > — 1 limit carefully. 



8.3.1 Can E' > Compensate for M' > 0? 



As noted in |]T^ the introduction of matter into a wormhole slows the progress of light rays 
through it. Can this effect be compensated for by a suitable choice of E' > 0? Since the 
vacuum case is M' = = E', for which we know the behaviour, we are only interested in the 
effects of varying M and E. 



We start with the gradient of the null rays, eq (|157| ) and use the substitutions ( pTI]) and 
(|113| ) with ( p.l4| ) for R, &l R' in terms of r] &l r, but we note that, if we choose the future AH 
— i.e. outgoing rays in a collapsing region — then 



j = +1 and Tx < 7] < 211. 



(191) 
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So the gradient of the null rays in terms of rj &i r \s 



dt 

dr 




M(l - COST]) 



smrj 



+ 1 



-/)Vi + /) 

3 sin 7/(77 — sin 77) \ /' 
2(l-cosr/)2 J 7 

sin r]{ri — sin rj) \ M' 



(1 — cosr/)^ 



E'' 

'e 



(1 - cosr/)2 y M 

Consider a region in which /?' > and M' > 0. Now since, in the above range of 77 

sinr/ 



(192) 



2 > 1 - COST/ > 



sin 77(77 — sin 77) 

1< 1 < OO 



> 



(1 — COS 77) 



> — CX), 



(1 — COS 77)2 



3 sin 77(7/ — sin 77) 
^ - ^ ^71 - ^' 

2(1 — COS7/j^ 



(193) 



the coefficient of M'/M is always positive, and the coefficient of E' / E is always negative. In 
particular, because of the no shell crossing condition (|97D, \E' / E\ < M'/3M, the E'/E term 
gives at most a partial cancellation of the M'/M term. Thus it is evident that varying E'/E 
cannot compensate for the effect of non-zero M' on the gradient of the radial rays. 

We turn to the AH equation (^) with (^) and(P^). We find the future AH equation 
in terms of i] and r may be written 



1 + 



(-/) 



sin?7 



\il + f) (l-COST/) 

3 sin 77(77 — sin 77) 



-ff/'a' 



smr] 



-f) 



2(1 — cos 77)2 
sin 77 



1 

M \ 



M \{l-cosr]Y^ 
sin ri{ri — sin 77) ^ 



(1 — cos 77)2 



\ (1 + /) (1-COS77) 



E' 

'e' 



(194) 



The solution is the parametric locus 7/ = rfAui'^)- 't is evident that if E' = 0, varying M' has 
no effect at all on the AH locus for a given M, as the solution is 



;-/) 



sm T] 



\{l + f) (1- cos 77) 



-1 



cos 77 = 1 + 2/ 



R = 2M. 



(195) 



Similarly E' has no effect when r] = n, and (see appendix [D]) when 77 = or 27r. On the other 
hand, the effect of varying E' is influenced by the value of M'. Analysing the slope of this 
curve leads to pretty daunting expressions, but is fortunately not necessary. 



Consider now the slope of a surface R{t,r) = aM{r) in a collapsing region 



R = aM 



R' 



2M 
oM 



+ f 



dt 
dr 



R! - aM' 



R=aM 



+ f 



(196) 
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where a > 0. This is null or outgoing timelike wherever 



{R' - aME'/E) R' - aM' 



For M' = 0, which forces E' = (by eq. (|97|)), the equality obviously requires a = 2, giving 
the event horizon in a vacuum model, and all R = aM surfaces are outgoing timelike for a > 2, 
and spacelike for a < 2. For M' > 0, the condition E' / E < M'/3M ensures the numerator 
of the Ihs is no less than 

R' - aM'/3. (198) 

For any given M' > 0, R' > and / > —1, this is always greater than the numerator on the 
rhs, so, to satisfy the equality, the denominator on the Ihs must be greater than that on the 
rhs, so once again 

a>2. (199) 

Thus R = aM surfaces can only be tangent to outgoing null rays for R = aM > 2M, and for 
R < 2M they are spacelike, incoming null, or incoming timelike. 

This allows us to conclude that, along the entire length of the future R = 2M surface, 
outgoing rays pass inside it, or run along it where M' = 0. By (0), the maximum R along 
any given constant r worldline in an elliptic region is at ?7 = tt, when 

Rmax = > 2M, (200) 

while R grows without bound in parabolic and hyperbolic regions. Thus every particle worldline 
encounters the future R = 2M surface (and the past surface), leaving no room for any rays 
that arrive at the future R = 2M surface to escape to J'^ . 

The time reverse of these arguments applies to the past AH, which lies in an expanding 
region and has incoming rays running along it or passing out of it. 

To complete the argument, we must consider the limits at the neck, / — ^ — 1 where several 
derivatives are zero. 



8.3.2 The AH at the Neck 

Now we turn to consider the AH at the neck. 

The differential of ( p.61| ) with ( |163| ) and dp = = dq gives us 



Q _ J R RE'\\Rdt + R'dr f dr 



v/TTT V E )\ R 2(1 + /) 

(R'dt + R"dr - - RMl^ + 

+R'dt + R"dr, (201) 
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dt 

dr 



AH 



+j 



R' 
R 



R! - 



RE'' 



vT+7 

R 



R" 



E 
R'E' 



Rf 



2(l + /)3/2 

RE" R{E'f 



R' - 



RE'" 



R' 



E E 
RE'' 



E^ 



E 



+ J 



R 



R' 



E 

_ RE' \ 



+ R' 



(202) 



At the neck of the wormhole, r = r„, the regularity conditions of sec. |^ give us the 
following limits, where Lf etc are being defined in each case: 

/ ^ -1, (203) 

/Vyi + Z ^ Lf' = V"/Lr > {Lf,f = 2/", (204) 

E'/^l + / -> Le' = 2E"/Lf. > 0, (205) 

R!l^lV] ^ L^' = 2R!'/Lf, > 0, (206) 

klJWf Lj,, = 2R"/Lf, = {d/dt)Ln'. (207) 



E' 0, 
R' ^ 0, 
i?' ^ 0, 

Thus the terms R'R'/^TTJ , R'RE'/E^/TTf , R'E'/E^/TTJ , R{E'f/E^^/TT7 and 
R' go to zero and the remaining numerator terms involving R cancel, down to 



dt 

dr 



-jR" 



AH,1<I 



dt 

dr 



J- 



2{R" -RE"/E) 



RE"' 



R 



VT+7 



R' 



RE'' 



(208) 
(209) 



Since the AH only intersects the neck when i? = 0, the behaviour of Rj \J\ + / and 
R? / y/1 + / must still be determined (and that of R! will be verified). 



At the moment of maximum expansion in the neck we have 

R = 2M , i? = , 



7] = 71 



COS?7 



and so, using 

R' = 



M 



;-/) 



—a 



'M' 



r 

sin?7 



(1 — cos rj) — 



'M' 



l + / = 0, 
sin rj = 0, 



3/'^ smri{ri — sinrj) 
2/, 



(210) 
(211) 



R" 



(1 — COST]) 

-1 M 



(1 - cos?7)2 (-/) 
sin rj ' 



[j] — sin r]) 
M 



'M' _ Sf 



+ 



(1 — cosri) 



(212) 



M 



7 



(1 — COS?]) 



(?7 — sin rf) 



-/) 



2/' fM' 9/' 
TIm ~ 8/ 



'M" 3/'' 



M { M" f" 2f' (M' f 



(213) 
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gives, by virtue of ( |204|) and ( |140|) : 

Lr, = 2{Lm' + MLf,) = 4(M" + Mf")/Lf, , B!' = 2{M" + Mf"). (214) 
We find the limit of R/ Vl + / at this point by combining (|161| ) and (|163| ), to obtain 



R 



-jL 



R' 



-j{M" + Mf'' 



MEN 



Lr, - RLe' /E M" + Mf" - ME" /E ' 



(215) 



VT+f 

so it is clear that W'E' j + / = 0. To check the limit of R! , the r derivative of (^) gives 

Rl 



1 (M' MR! , f\ RR' [Lm' MLr, , Lf, 



R\R R^ 



+ 



R 



i?2 + 2 



(216) 



and because of (|214| ) and R = 2M all terms in the bracket cancel, verifying that RR! / ^/l + / 
0. 



These together with (^ give us 

dt 
dr 



yl_H',MEN 



dt 
dr 



n,MEN 



AMf"{M" + Mf") 

' Lf, {M" + Mf" - ME" IE) ' 
A{M" + Mf" -ME" IE) 



(217) 
(218) 



For a light ray to pass through the neck at the moment of maximum expansion without 
falling inside the AH, we need (it/dr |^ji/^men > '^V'^'^U.men. in other words 



dtldr 



yl_H',MEN 



Mf"{M" + Mf" 



dtldr\n,M^^ (M" + Mf" - ME"IE)^ 



>1, 



(219) 



MEN 



or 



(M" + Mf") - JMf"{M" + Mf") E" {M" + Mf") + JMf"{M" + Mf") 

< — < ^ . 

MEM 

(220) 

Since M(r) is positive, and both M(r) and /(r) are at a minimum at the neck, i.e. M" > 



&i f" > 0, we have ^Mf"{M" + Mf") < M" + Mf", and so both upper and lower limits 
are real and positive. 

Can this requirement be satisfied without creating shell crossings? The only relevant 
condition is the one for e = +1, i?' = 0, / = -1, R" > 



E" 



< 



M" 
3M' 



(221) 



To be able to satisfy this as well as ( |220|) we would need 



M" + Mf" - \lMf"{M" + Mf") M" 



M 



3M' 



(222) 
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but this leads to 

M"(4M" + 3Mf) < 0, (223) 

which is clearly not possible. Indeed, although p is zero rather than divergent where E' / E = 
M'/3M, where E'/E exceeds M'/3M, the density is negative at all times. 

Putting the maximum value, E"/E = M" /3M into ( PD gives 



dt/dr\ 



AHMEN 



dt/ dr\n,uET<s 



9MnM" + Mn (224) 



which rises from at f" /M" = 0, and asymptotically approaches 1 as f"/M" oo, i.e. 
vacuum. 

Therefore, even at the neck, E' > cannot compensate for M' > 0, and all rays passing 
through this event remain within R < 2M, passing from inside AH^ to inside AH+. 



8.3.3 Summary 

In a Szekeres wormhole, every particle worldline encounters R = 2M, twice for most r values 
and once where / = —1, making this a pair of 3-surfaces that span the spacetime. The 
apparent horizons coincide with R = 2M at an extremum of R{t = const, r) — a neck or belly 
— where / = —1. Where M' = (vacuum), the R = 2M surfaces are (locally) null. 

Assuming there is matter (M' > 0) somewhere within the elliptic region describing the 
neck, and assuming the two regions, r ±oo, are asymptotically flat, i.e. M Mtot = con- 
stant {E' 0), then the event horizon is the set of rays that are asymptotic to i? = 2M, 
but always lie outside. The future event horizon EH+ emerges from R = 2M surface, and 
vice- versa for EH~. 

Thus we conclude that the causal structure of a regular Szekeres wormhole is only a quan- 
titative modification of the LT wormhole (dense black hole), and the possible causal diagrams 
for Szekeres models are essentially the same as those for LT models, as given in ||T^. 



8.3.4 Numerical Examples 

A few numerical examples were produced as follows. 

We choose the 3 LT arbitrary functions to produce a Kruskal-like topology, with the neck 
at r = 0, that is mirror symmetric about r = and t = 0. The choice must therefore satisfy 

/(O) = -1, /'(O) = 0, /" (0) > 0, M'(0) = 0, M"(0) > 0, a(r) = -6(r); 

M = Mo{l + Miry , Mo,Mi> 0, (225) 
/ = -exp(-rVr.) , r, > 0, (226) 
a = -TTM/{-ff\ (227) 

We want to choose E to maximise the effect of E' ^ along one particular radial path. By 
setting 

P(r) = = Q(r), (228) 
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so that (0) is 

the maximum E' / E becomes 

(230) 



E' 




S' 




max 


s 



along the direction (p, g) = (0,0), i.e. 6' = 0. Since numerical integrations will only be done 
along this path and the 9 = n one, we treat E as a function of r only. We make E' / E as large 
as possible without violating the no shell crossings condition E' / E < M' /3M with 

E = Eo{l + Eir^) + E2 , Eo,Ei,E2>0, (231) 

where the shell crossing occurs somewhere if E2 = 0. 

Because of the two reflection symmetries, we can start integrating a null ray from maximum 
expansion at the neck, 

r] = 7r , r = 0, (232) 

where AH+ &l AH~ meet. The symmetry means that integrating forwards along increasing r 
&L t and integrating backwards along decreasing r &l t \s the same thing, so one integration 
actually traces both halves of the same ray. Rays that don't pass through this point require 
two separate parts to the integration, one from maximum expansion towards r &l t increasing, 
and the other towards r &l t decreasing, with careful treatment of the neck limits where r goes 
through zero. 

The following runs were done: 

• Test 1 — the vacuum case: 

Mo = 1 , Ml = , ^0 = 1, ^1 = 0, E2 = , r, = 1. (233) 

As expected, we found that the fast AH, the slow AH, the fast null ray, and the slow null 
ray were all the same. 

• Test 2 — the LT case: 

Mo = l, Ml = 0.1, Eo = l, Ei=0, ^2 = 0.01, r, = 1. (234) 

Here the fast &l slow rays were the same, and the fast &l slow AHs were the same, but 
the rays fell inside the AHs, as expected. 

• Test 3 — a Szekeres version of above LT case: 

Mo = 1 , Ml = 0.1 , ^0 = 1, ^1 = 0.1 , E2 = 0.01 , r, = 1. (235) 
The AHs and rays were split on either side of the Test 2 curves. 

• Run 1 — medium M'/f 

Mo = 1 , Ml = 1 , Eo = l, ^1 = 1, E2 = 0, r, = 1. (236) 
We found that the rays &l AHs were well split, while the rays were strongly trapped. 
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• Run 2 — low M'lf 

Mo = l, Ml = 1 , ^0 = 1, Ei = l, ^2 = 0.1, r, = 0.01. (237) 

The rays were mildly split, the AHs were indistinguishable in the range plotted, and the 
rays were mildly trapped. 

• Run 3 — slightly less low M' / f 

Mo = 1 , Ml = 2 , ^0 = 1, ^1 = 2, E2 = 0.01 , r, = 0.01. (238) 
This was very similar to the previous run, with the rays less mildly trapped. 

• Run 4 — high M' / f 

Mo = 1 , Ml = 3 , ^0 = 1, ^1 = 3, E2 = 0.0001 , r, = 10. (239) 
Here the rays &l AHs were well split, and the rays were very strongly trapped. 

These examples cover the main possibilites, and run 1 is shown in Fig. ^. 

9 Conclusions 

Szekeres (S) models are a generalisation of the spherically symmetric LemaTtre-Tolman (LT) 
models. Both describe inhomogeneous dust distributions, but the former have no Killing vectors. 
There are 3 arbitrary functions of coordinate radius in LT models (M, / Su. a), and a further 3 
in S models [S, P &i Q). 

For quasi-spherical Szekeres (S) models, we established 3 sets of regularity conditions — 
the conditions for a regular origin, the conditions for no shell crossings, and the conditions for 
regular maxima and minima in the spatial sections. The last two contain exactly those for the 
LT models, with extra conditions on the arbitrary functions that are peculiar to S. Thus, for 
every regular LT model that is non-vacuum (M' > 0) at least somewhere, one can find regular 
S models that are anisotropic versions of the same topology. (For vacuum, M' = 0, S models 
must be spherically symmetric.) 

Since LT models can reproduce the Schwarzschild-Kruskal-Szekeres topology of a wormhole 
connecting two universes, but with non-zero density everywhere, this is also possible with 
S models. In the vacuum case (M' = 0) this gives the full Kruskal manifold in geodesic 
coordinates. It is known that the presence of matter in such models inhibits communication 
through the wormhole and splits the event horizons. We investigated the S wormhole models, 
considering apparent horizons and the paths of 'radial' null rays, which, while not geodesic, are 
the fastest paths out of a wormhole. We showed that, even though the S model's anisotropy 
makes the proper separation of consecutive shells shorter along certain directions, and null 
motion faster along those same directions, this is not enough to compensate for the retarding 
effect of matter. Thus the causal structure of an S wormhole is the same as that of the 
corresponding LT model. 

We also considered whether the two universes on either side of a wormhole could be joined 
across a 3-surface, making a handle topology. It was found that a smooth junction is not 
possible at any finite distance, as a surface layer would be created. This conclusion applies to 
LT models and to the vacuum case — a Schwarzschild wormhole — too. 
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Fi gure 6; The (r — t) diagram for the Szekcrcs model defined for run 1, showing the fast 
and slow future apparent horizons (fA+ and sA+), and past apparent horizons, the fast 
and slow rays that pass through O — the neck at the moment of maximum expansion — 
towards r increasing (fR+ and sR+), and rays through O going towards r decreasing, as 
well as rays going through other points. T is the moment of time symmetry which is also 
the simultaneous time of maximum expansion, and N is the locus of the neck r = 0. Note 
that fA+ & sA+ arc two different intersections of the future apparent horizon AH^ in two 
different radial directions — the fast & slow poles where E' / E takes extreme values. Note 
also that there is no origin R{r = ro,t) = in wormhole models. 
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A The hypersurface of zero mass-dipole intersects every 

{t = const, r = const) sphere. 

This hypersurface is given by 

E'/E = {A' + C')/{A + C). (240) 

Since E'/E at constant r is bounded (see (|64D), it must be verified whether eq. ( |240|) has a 
solution in every sphere of constant t and r. The solution will exist when 

(E'/EU^ < (A' + C')/iA + C)< (^7^)max. (241) 
Since {E' / E)^^^, = -(-E7-E)max, eq. (^) is equivalent to 

{A' + C'y/{A + Cf < (i?7^)extrcme = ^2 + + • (242) 

We have 

A + C = ^(l + + + S^) 

A' + C' = ^[S^-P^-Q''-l) + ^ {PP' + QQ') . (243) 
Substituted in ( |242| ), this leads to 

AS^S'^ (l + + - 455' {PP' + QQ') (s^ -P^-Q^-l) 

-4S\PP' + QQ'f + (P'^ + g'^) (l + P2 + g2 + S^) > 0. (244) 
The discriminant of this with respect to S' is 



2 r 



A = -16S^ (l + p2 + g2 + {PQ' - QP'f + + g'2 ^ (245) 
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and is always negative unless P' = Q' = 0. This means that with {P',Q') ^ (0,0), the l.h.s. 
of ( p44| ) is strictly positive. Even when P' = Q' = 0, it is still strictly positive unless S" = as 
well. However, P' = Q' = S' = implies A' = C = Q and £" = on the whole sphere, and 
then the dipole component of density Ap = 0; i.e. on such a sphere the density is spherically 
symmetric. Hence, apart from the spherically symmetric subcase, eq. (|241| ) is fulfilled, with 
sharp inequalities in both places. This means that the Ap = hypersurface intersects every 
{t = const, r = const) sphere along a circle parallel to the E' = circle (see remark after eq. 

my 



B Matching the Szekeres Metric to Itself 

We here lay out the calculations necessary for matching the Szekeres metric across a comoving 
surface to some other metric, and in particular to another Szekeres metric. 

Given a comoving surface, 

rs = Z{p, q) (246) 

and surface coordinates, 

e = (t,p,g) (247) 
we calculate the basis vectors in the surface, 

= (248) 



the 1st fundamental form, 



the normal vector. 



, n^n'' = 1 , n^ef = 0, (250) 

and the 2nd fundamental form 



-nt (- 



d'^x^ „\ dx^ dx 



Using GRTensor/GRJunction pD| we find the following for the intrinsic metric: 



hi = -1, (252) 

(253) 



Z^{R' - RE'/EfE^ + R\e + /) 
yw - E\e + f) 



ZpZg{R' - RE'/Ef 

^TTO ' ^'''^ 



Zl{R-RE'/EfE^ + R\e + f) 
= E^^^ , (255) 
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the surface basis vectors: 



the surface normal: 



rir — 



where 



el = (1,0,0), 

e? = (0,1,0), 
ef = (0,0,1), 



{R'E - RE')R 



EA 

Zpi^R'E — RE'^R 

Zq{R'E - RE')R 

'eK ' 

Z 



Z 



A = (R\e + f) + {Zl + Zl){IiE-RE'ff\ 



and the extrinsic curvature: 

[Zp{RR' - RR')] 



pt 



[Zq{RR' - RR')] 



(256) 
(257) 
(258) 
(259) 



(260) 
(261) 
(262) 
(263) 
(264) 

(265) 
(266) 



K — 
pp — 



\ 2RE{R'E - RE'){e + f)Z^ 



pp 



2E\e-rf)A\ 

+2(R'E - RE'f{EpE' - E'pE)Z^ 
+2(R'E - RE'f{EqE' - E'gE)Z^Zq 

+ {2 [SER'RE' + RR"E'^ - {E'fR^ - 2E^{R'f - R^E"E 

- Rf'E{R!E - RE')]zl 
-2R{2RE'^E - EpR'E - REpE'){e + f)Zp 
-2REq{R'E - RE'){e + f)Zg 

-2R\e + fA, 



(e + f) 



(267) 



K. 



pq 



\ 2RE{R!E - RE'){e + f)Z^ 



pq 



2E^{e + f)A\ 
+2{R'E - RE'f{EpE' - E'pE)Z^Zg 
+2{R'E ~ RE'f{EqE' - E'qE)ZpZl 

+ {2 ^ER!RE' + RR"E'^ - {E'fR^ - 2E^{R'f - R^E"E 
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- RfE{R'E - RE')}ZpZq 
-2RE{RE'^-R'Eg){e + f)Zp 

-2RE{RE'p-R'Ep){e + f)zJj, 



(268) 



+2{R'E - RE'f{EpE' - E'^E)ZpZl 
+2{R'E - RE')\EgE' - E'^E)Z^ 

+ {2 [SER'RE' + RR"E^ - {E'fR^ - 2E^{R!f - R^E"e] (e + /) 

- Rf'E{R'E - RE')}ZI 
-2REp{R'E - RE'){e + f)Zp 
-2R{2RE'gE - EgR'E - REqE'){e + f)Zq 

-2R'{e + fA, (269) 



where all quantities are evaluated on E. 



C The Acceleration of a Given Tangent Vector 



Starting from 

the individual acceleration components for a given /c" in the Szekeres metric are 

a* = k^dpk' + r'rr{k')^ + r\p{kP)^ + r\g{k'if 

= k^dtk^ + k'^drk^ + k^dpk^ + kWgk^ 



(270) 
(271) 



+ {R 



, RE'\(^, RE'\ 1 



E 



R'- 



E / e + f 



RIV 
E^ . 



{k^f + {k'^f) , (272) 



a — 



+rpp(F)' + r,,(A;«)2 



k'dtk' + k'drk' + k^dpk' + k'^dgk' + 



k'k 



r 



(R' - if) 2(6 + /) 



2R _ -^'-^p 
(D^ - (^f _ J) ^ fe^^^ 



^2 (i?' - ^) 



(273) 



(274) 
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oF = A;^9/9F + 2r^'tpA;*F + r^,(A;")2 + 2r%A;''F 



= A;*atF + k'drkP + k^dpkP + k'^dgk^ + ^k^k^ 



2R, 
~R 



i?(e + /) 



~l~r^pp(/c^) + 2V^pqk^k'^ + r^q,q(A:^) 



2 /? 

A;*9tA;« + A;''9^A;« + F6>pA;« + k'^dqk'^ + —k^k'^ 

R 



R 



R(e + f) 

For "radial" paths k^ — — k"^, dpk"' = = dqk"' these reduce to 

a* = k'dtk' + k''drk^ 



E 



E J (e + f) 



2 R'- 



I RE' 



k'k' 



+ 



{R' - ^) 

(i?/ _ if) 2(6 + /) 

WW) 



R{e + f) 



Using the "radial" null condition 



k 



t _ jk 



' 'r' - ^ 
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the acceleration becomes 



J 



J 



R'- 



RE'^ 



E 



R' 



k'k' + R' - 



f 



RE'' 



+ R'- 



RE'' 



RE'V 
e) 2(e + /) 

RE'' 



R' 



RE'' 



k'-drk' } + {R' - 



E 



R'- 



E J 
RE' 



1 



J 



+ 



R 



I RE' 



R' 



RE'' 



¥dtk'' + 2 i?' 



E J (e + f) 
RE' 



(kr, (284) 



r\2 



E 



(kl 



r\2 



(kl 



{R' - if) 2(e + /) 



R'-^) {EE'^-E'E,) ^^^^^^ 



{k^f. 



R' - ^) {EE'^ - E'E,) 
R{e + f) 



(285) 

(286) 
(287) 



While E' = gives the expected LT values, we note that E{r,p, q) determines whether and 
a'^ are zero or not. 

In the e = +1 case, by (|16|), ( |60D and (|62|) , the extremes of on a given 2-sphere are 
located at 



Pe 



P + 

Q + 



P'S 



( ± ^{S')^ + {P'f + {Q'f ~S') 

Q^s 

( ± ^{S'Y + {P'Y + {Q'Y - s' 



(288) 
(289) 



and it is easily verified that = = a'^ in these two antipodal directions. It follows that 
initially radial geodesies in these directions remain radial if &l are constant with r. For 
example, if Pe = = qe, this would require arbitrary functions satisfying 



P' 



2S'S 



Q' 



or 

or 
or 



P S^-P^-Q^ Q 
P = . «' '''' 
P' 2S'S 



when P' = 



p g,_p, , Q = , when Q' = 
P = , (5 = 0, when P' = and Q' = 



(290) 

(291) 

(292) 
(293) 
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D Other Features of the AH 



The FLRW Case 

The dust FLRW limit is M = Mor^ / = -kr^, a = 0, R = rS{t), and the R = 2M locus is 
given by 

SitAH) = (294) 
In the collapse phase of a A; = +1 model, the time of the future AH is 

COST] AH = l-2r^ (295) 
tAH = Mo[tt + arccos(2r2 - 1) + 2rVl - ] (296) 

which has slope 



dt\ 4MoH 



r 



2 



while the light rays have slopes 

fdt\ S 2Mor^ 



(297) 



(298) 



Clearly the future AH is incoming timelike. The result extends to all k values, and the converse 
holds for the past AH in the expansion phase. 



Behaviour Near the Bang or Crunch 



Consider eq (|194| ) for the locus of the future AH in a collapsing elliptic region, f < 0, tt < r] < 
27r, in terms of parameter i]. Near the crunch, rj = 2n — r] ^ 0, we find 

r / 127r\ /87r\ 1 



+ 




\ (1 + /) 



(299) 



As noted previously, when E' = 0, the solution makes the first bracket zero, fj ^ 2y/—f 
(Even in this case, where we know R = 2M is the AH, the fact that R' diverges at i? = 
means we must multiply through by r] to make the rhs zero there.) Notice too that the no shell 
crossing condition ( [126| ) ensures the second bracket is generically non-negative where M' > 
and non-positive where M' < 0. Assuming we aren't near an origin, < M < oo, it is clear 
that, even if / 0, the last two terms in this second bracket are divergent, with the middle 
one dominant, making the E' term negligible. Thus rj ^ 2-^/— / is still the solution in the 
limit. However (|T2|) in the rj 
to 



limit shows that the time from AH to crunch goes 



b - tAH 



M 



AM 



(300) 



(-/)3/2 6 3 

where the crunch time b{r) is defined in ( |128D . Therefore the future AH does not intersect the 
crunch away from an origin. The result is just the time reverse for the past AH near the bang 
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in an expanding region, and a similar calculation applies for hyperbolic or extended parabolic 
regions, giving the same result. 

An / = locus is where an interior elliptic region joins to an exterior hyperbolic or parabolic 
region. The transition involves the lifetime of the worldlines diverging, so either the crunch 
goes to the infinite future, or the bang goes to the infinite past. Since there is only one AH 
in a hyperbolic or parabolic region, one of the two AH loci in the elliptic region also exits to 
infinity before / = is reached. 

A third possibility where / is only asymptotically zero (the asymptotically flat case) is that 
both the bang and the crunch diverge to the infinite past and future, and the two AHs go with 
them. 



Behaviour Near an Origin 



Consider eq ( p.94|) again. Near a regular origin, along a constant t or constant i] surface (see 
section |), M ^ /)^/^, E ^ /)"^^ for some positive constants ^ &l v, and / ^ 0, so 
that 







1 + V-/ 

/ 



(1 — COST]) 

sin 7] 



smrj 



jj, \{1 — cosriy 



+ 



2/ 



(1 — cosrf) 
We divide through by /' and define X = 



(301) 



-/ which is positive for r] > n, giving 







SUIT] 



X n 



/if V(l-cos77)V 2/ 1-X2/ 



(302) 



Though a' and M' always have opposite signs, /' may have either sign in an elliptic region, but 
in general we don't expect terms to cancel in the curly brackets. Thus, whether or not the a' 
term diverges, we must have X ^ \, i.e. sin// so that 

= 27r - r/ ^ 2^f. (303) 
Unlike the previous case, though, M ^ ensures the AH intersects the crunch here. 



h - tAH /i 



0. 



(304) 



As always, the time reverse applies in an expanding phase, and the hyperbolic and parabolic 
cases give the same result. 
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